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Abstract 

This work gathers new results concerning the semi-geostrophic 
equations: existence and stability of measure valued solutions, exis- 
tence and uniqueness of solutions under certain continuity conditions 
for the density, convergence to the incompressible Euler equations. 
Meanwhile, a general technique to prove uniqueness of sufficiently 
smooth solutions to non-linear ly coupled system is introduced, using 
optimal transportation. 
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1 Introduction 

The semi-geostrophic equations are an approximation to the Euler equa- 
tions of fluid mechanics, used in meteorology to describe atmospheric flows, 
in particular they are believed (see jT^j) to be an efficient model to describe 
frontogenesis. Different versions (incompressible P , shallow water ^U] , com- 
pressible [TT]) of this model have been studied, and we will focus here on the 
incompressible 2-d and 3-d version. The 3-d model describes the behavior 
of an incompressible fluid in a domain O C I 3 . To the evolution in Q is 
associated a motion in a 'dual' space, described by the following non-linear 
transport equation: 

d t p + V • (pv) = 0, 
v = (W(x) -x) ± , 
det D 2 ^ = p, 
Pit = 0) = p°. 

Here p° is a probability measure on R 3 , and for every v = (vi,V2,v 3 ) £ R 3 , 
V" 1 stands for (— v 2 , v 1, 0). The velocity field is given at each time by solving 
a Monge- Ampere equation in the sense of the polar factorization of maps (see 
jH]), i.e. in the sense that \& is convex from R 3 to R and satisfies V\l/#p = 
Xn£ 3 , where £ 3 is the Lebesgue measure of R 3 , and xn is the indicator 
function of Q. For compatibility Q has Lebesgue-measure one. This model 
arises as an approximation to the primitive equations of meteorology, and 
we shall give a brief idea of the derivation of the model, although the reader 
interested in more details should refer to |12j . 

In this work we will deal with various questions related to the semi- 
geostrophic (hereafter SG) system: existence and stability of measure-valued 
solutions, existence and uniqueness of smooth solutions, and finally conver- 
gence towards the incompressible Euler equations in 2-d. The results are 
stated in more details in section 1. 
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1.1 Derivation of the semi-geostrophic equations 

We now give for sake of completeness a brief and simplified idea of the deriva- 
tion of the system, inspired from [lj, and more complete arguments can be 
found in IT2l. 



Lagrangian formulation 

We start from the 3-d incompressible Euler equations with constant Coriolis 
parameter / in a domain Q. 

Tj7 + / v± = -Vp - Vy>, 
Dt p 

Dp 

V.v = 0, ^ = 0, 
v ■ dQ = 0, 

where — stands for d t + v • V, and we still use = (— t> 2 ,t>i,0). The 

term V<£> denotes the gravitational effects (here we will take ip = gx% with 

constant g), and the term fv x is the Coriolis force due to rotation of the 

Earth. For large scale atmospheric flows, the Coriolis force /v -1 dominates 
Dv 

the advection term —— , and renders the flow mostly bi-dimensional. We use 

Dt' J 
the hydrostatic approximation: d X3 p = —pg and restrict ourselves to the case 

p = l. 

Keeping only the leading order terms leads to the geostrophic balance 

that defines v g , the geostrophic wind. Decomposing v = v 9 + \ ag where 
the second component is the ageostrophic wind, supposedly small departures 
from the geostrophic balance, the semi-geostrophic system reads: 

V • v = 0, 

where V h = (9 Xl , d X2 , 0). Note however that the advection operator d t + v- V 
still uses the full velocity v. Introducing the potential 

$ = \\*h\ 2 + r 2 v, 

with iff = (xi,a;2,0), we obtain the following 

^V$(t,x) =/(x-V$(t,x)) ± . 
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We introduce the lagrangian map g : x M + i— > Q giving the position at time 
t of the parcel located at Xq at time 0. The previous equation means that, 
if for fixed iGflwe consider the trajectory in the 'dual' space, defined by 
X(t,x) = V$(t,g(t, x)), we have 

d t X(t,x) = f(g{t,x)-X(t,x)) ± . 

By rescaling the time, we can set / = 1. As stated the system looks under- 
determined: indeed $ is unknown; however we have the condition X(t,x) = 
V<E>(i, g(t, x)). Moreover, the dynamic in the x space being incompressible 
and contained in Q, the map g(t, •) must be measure preserving in Q for each 
t, i.e. 

£*(g(ty 1 (B)) = £\B) 

for each B C £1 measurable (where C 3 denotes the Lebesgue measure of M 3 ). 
We shall hereafter denote by G(Q) the set of all such measure preserving 
maps. Then Cullen's stability criteria asserts that the potential $ should be 
convex for the system to be stable to small displacements of particles in the 
x space. Hence, for each t, <3> must be a convex function such that 

x(t,.) = v$(t,g(t,.)), 

with g(t, •) G G(Q). In the next paragraph we shall see that, under very 
mild assumptions on X, this decomposition, called polar factorization, can 
only happen for a unique choice of g and V$. Now if $* is the Legendre 
transform of $, 

$*(?/) = sups • y — 

then V$ and V$* are inverse maps of each other, and the semi-geostrophic 
system then reads 

D X 

— = (V<T(X(t))-X(t))\ 
V$*{t)oX(t) e G(Q). 

In the next paragraph, we expose the results concerning the existence and 
uniqueness of the gradients V$, V$*. 

1.2 Polar factorization of vector valued maps 

The polar factorization of maps has been discovered by Brenier in [3] . It has 
later been extended to the case of general Riemannian manifolds by McCann 
in j20!. 
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The Euclidean case 

Let Q be a fixed bounded domain of M d of Lebesgue measure 1 and satisfying 
the condition C d (dVl) = 0. We consider a mapping X e L 2 (fi;M d ). We will 
also consider the push- forward of the Lebesgue measure of fl by X, that we 
will denote by X#xn£ d — dp (or, in short, X#dx) and which is defined by 



Let V be the set of probability measures M d , and V 2 the subset of V where 
the subscript a means absolutely continuous with respect to the Lebesgue 
measure (or equivalently that have a density in L 1 (M d )), and the superscript 
2 means with finite second moment, (i.e. such that 



Note that for X e L 2 (Q,M. d ), the measure p = X#dx has necessarily finite 
second moment, and thus belongs to V 2 . 

Theorem 1.1 (Brenier, |3J). Let Vt be as above, X 6 L 2 (tt;M. d ) and p = 
X#dx. 

1. There exists a unique up to a constant convex function, that will be 
denoted $[p], such that: 



2. Let ^[p] be the Legendre transform o/$[p], if p € V 2 , ^[p] is the unique 
up to a constant convex function satisfying 







3. If p (zV 2 , X admits the following unique polar factorization: 



X = V$[p]og, 



with $[p] convex, g measure preserving in Q. 
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Remark: \Ef [/>],<& [jo] depend only on p, and are solutions (in some weak 
sense) respectively in M. d and Q, of the Monge- Ampere equations 



When \1/ and $ are not in Cf oc these equations can be understood in the 
viscosity (or Alexandrov) sense or in the sense of Theorem ILlj which is 
strictly weaker. For the regularity of those solutions and the consistency of 
the different weak formulations the reader can refer to jHj. 

The periodic case 

The polar factorization theorem has been extended to Riemannian manifolds 
in [20] (see also for the case of the flat torus). In this case, we consider 
a mapping X : M. d i— »■ M. d such that for all p G Z d , X(- + p) = X + p. Then 
p = X#dx is a probability measure on T d . We define *&[p), $[p] through the 
following: 

Theorem 1.2. Let X : M. d — > M. d be as above, with p = X#dx. 

1. Up to an additive constant there exists a unique convex function 
such that $[p](x) — x 2 /2 is Z d -periodic (and thus V$[p](x) — x islP 
periodic), and 



2. Let ty[p] be the Legendre transform o/$[p]. If p is Lebesgue integrable, 
^[p] is the unique up to a constant convex function satisfying 
^![p\(x) — x 2 /2 is Z d -periodic (and thus V^[p](x) — x is Z d periodic), 
and 



3. If p is Lebesgue integrable, X admits the following unique polar factor- 



det D 2 ^ = p, 
p(V$) detD 2 $ = 1. 





ization: 



X = V$[p] o g 



with g measure preserving from T d into itself, and $[p] convex, $[p] — 
\x\ 2 /2 periodic. 
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Remark 1: From the periodicity of V$[p](x) — x, V^[p](x) — x, for every 
/ Z d -periodic, /(W[p]), /(V$[p]) are well defined on R d /Z d . 

Remark 2: Both in the periodic and non periodic case, the definitions 
of \I/[p] and $[p] make sense if p is absolutely continuous with respect to 
the Lebesgue measure. If not, the definition and uniqueness of $[p] is still 
valid, as well as the property V$[p]#p = Xn£ d - The definition of \l/[p] as 
the Legendre transform of $[p] is still valid also, but then the expression 

f(V^[p](x)) dp(x) does not necessarily make sense since V\l/ is not neces- 
sarily continuous. Moreover the polar factorization does not hold any more. 
Remark 3: We have (see |§|) the unconditional bound 

||W[p](:r) - xWloo^ < Vd/2 
that will be useful later on. 

1.3 Lagrangian formulation of the SG system 

From Theorems II. 1[ 11.21 the Lagrangian formulation of the semi-geostrophic 
equation then becomes 

— = [V*(X)-X]\ (1) 
tf = tt[p], p = X#dx. (2) 

1.4 Eulerian formulation in dual variables 

In both cases (periodic and non periodic) we thus investigate the following 
system that will be referred to as SG: we look for a time dependent proba- 
bility measure t — ► p(t, ■) satisfying 

d t p + V • (pv) = 0, (3) 
v(t,x) = (W[p(t)]Or)-a;)\ (4) 
p(t = 0) = p°. (5) 

Weak solutions (which are defined below) of this system with L p initial data 
for p > 1 have been found, see [l], [101, [IS] . 

1.5 Results 



In this work we deal with various mathematical problems related to this sys- 
tem: we first extend the notion of weak solutions that had been shown to exist 
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for p e L°°(R + ,L g (M 3 )), q > 1 ([[], [TO]), and then for p G L 1 ^ 3 )) 
(|18j), to the more general case of bounded measures. The question of ex- 
istence of measure-valued solutions was raised and left unanswered in those 
papers, and we show here existence of global solutions to the Cauchy prob- 
lem with initial data a bounded compactly supported measure, and show the 
weak stability/compactness of these weak measure solutions. 

Then we show existence of continuous solutions, more precisely, we show 
local existence of solutions with Dini-continuous (see (JEJ)) density. For this 
solutions, the velocity field is then C 1 and the Lagrangian system (I1I2|1 is 
defined everywhere. 

We also show uniqueness in the class of Holder continuous solutions (a 
sub-class of Dini continuous solutions). This proof uses in an original way 
the optimal transportation of measures by convex gradients and its regularity 
properties, and can be adapted to give a new proof of uniqueness for solutions 
of the 2-d Euler equation with bounded vorticity, but also for a broad class 
of non-linearly coupled system. The typical application is a density evolving 
through a transport equation where the velocity field depends on the gradient 
of a potential, the potential solving an elliptic equation with right hand side 
the density. Well known examples of such cases are the Vlasov-Poisson and 
Euler-Poisson systems. 

Finally, in the 2-d case, we study the convergence of the system to the 
Euler incompressible equations; this convergence is expected for p close to 
1, since formally expanding \1/ = x 2 /2 + eip, and linearizing the determinant 
around the identity matrix, we get 

det£> 2 ^ = l + eA^ + 0(e 2 ), 

and the Monge- Ampere equation turns into the Poisson equation 

e 

After a proper time scaling, p satisfies 

d t p + V- (^VV), 
&ii = p, 

that we recognize as the vorticity formulation of the 2-d Euler incompressible 
equation. The study of this 'quasi-neutral' limit is done by two different ways: 
One uses a modulated energy method similar as the one used in jl] and [3] 
and is valid for weak solutions. The other uses a more classical expansion 
of the solution, and regularity estimates, and is similar to the method used 
in [123. The second method yields also a time of existence for the smooth 
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solution that goes to infinity, as the scaling parameter e goes to 0. From a 
physical point of view, this asymptotic study may be seen as a justification 
of the consistency of the semi-geostrophic approximation. 

2 Measure valued solutions 

2.1 A new definition of weak solutions 

We have first the following classical weak formulation of equation (J3|): 
p G C(JRL+, L X (]R 3 ) — w) is said to be a weak solution of SG if 

VT > 0, G C™([0,T] x R 2 ), 
/ d t ipp + Vip- (W[p] - x)- 1 p dtdx 



where for all t, $ [p] is as in Theorem 11.11 The problematic part in the case 
of measure valued solutions is to give sense to the product pV\I/[p] since at 
the point where p is singular V^fp] is unlikely to be continuous. Therefore 
we use the Theorem II. II to write for any p G V 2 (M. 3 ) 



(the integrals would be performed over T 3 in the periodic case). The property 
V$[p]#Xq£ 3 = p is still valid when p is only a measure with finite second 
moment (see Remark 2 after Theorem II .2|) . Therefore, the formulation on 
the right hand side extends unambiguously to the case where p ^ L 1 (IR 2 ). 

Geometric interpretation 

This weak formulation allows has a natural geometric interpretation: at a 
point where \l/[p] is not differentiable, and thus where <9\I/[p] is not reduced 
to a single point, V^fp] should be replaced by d^/[p] the center of mass of 
the (convex) set c^p]. 

This motivates the following definition of weak measure solutions 

Definition 2.1. Let, for all t G [0, T], p(t) be a probability measure o/R 3 . 
It is said to be a weak measure solution to SG with initial data p° if 

1- The time dependent probability measure p belongs to C([0, T], V — w*), 
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2- there exists t — ► R(t) non-decreasing such that for all t G [0, T], p(t, ■) 
is supported in B(0,R(t)), 

3- for allT>0 and for all <p e C c °°([0, T] x M 3 ) we have 
d t (p(t,x) dp(dt,x) (6) 

+ / V<p(t, V$[p(t)](x)) • x 1 - dtdx - Vip(t,x) ■ x L dp(dt,x) 

J[0,T]xn J[0,T]xR3 

(p(T,x)dp(T,x) dx — I <p(0, x)dp°(x) dx. 



[0,T]xR 3 



This definition is consistent with the classical definition of weak solu- 
tions if for all t, pit, ■) is absolutely continuous with respect to the Lebesgue 
measure. 



2.2 Result 

Here we prove the following 

Theorem 2.2. 1. Let p° be a probability measure compactly supported. 
There exists a global weak measure solution to the system SG with 
initial data p° in the sense of Definition \2. li 

2. For any T > 0, if (p n ) nt =n is a sequence of weak measure solutions on 
[0, T] to SG with initial data (p°) ng N ; supported in B R for some R > 
independent of n, the sequence (p„)„ e N is precompact in C([0,T],V — 
w*) and every converging subsequence converges to a weak measure 
solution of SG. 

Proof of Theorem 12.21 

We first show the weak stability of the formulation of Definition (j2.1|) . and 
the compactness of weak measure solutions. We then use this result to ob- 
tain global existence of solutions to the Cauchy problem with initial data a 
bounded measure. 

Weak stability of solutions 

We consider a sequence (p„)„ e N of solutions of SG in the sense of Definition 
12.11 The sequence is uniformly compactly supported at time 0. We first show 
that there exists a non-decreasing function R(t) such that p n (t) is supported 
in B(R(t)) for all t,n: 
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Lemma 2.3. Let p G C([0, T], V(M. 3 ) - w*) satisfy ffl, let p° = p(t = 0) be 

supported in B(0, R°), then pit) is supported in B(0, R°+Ct) , C = sup yen {\y\}. 

Proof. Consider any function £ e (i, r) G C°°(R) such that 

6(0, r) = 1 if - oo< r < 
6(0, r) = if r > i2° + e, 
&(f,r)=&(r-Ct), 

with £(0, •) non increasing. Then compute 
— / £ e (t, M) dp(t,x) 

= - J d r Ut)x\)C d P (t,x) + J^d r Ut)V^ P (t)}\)^^-X ± dx 



> / d r Ut,\V$[p(t))\)(-C+\x\) dx 
Ji 

> 



n 



since, by definition, for x G f2, Ixl < C and £ is non increasing with respect 
to , Note a ls0 that we have Ud"/ V.R( t , W) ] • * S 0. We 

know on the other hand that 

&(0,|x|)dp(0,x) = l, 
&(t,|z|)dp(t,z) < 1, 

therefore we conclude that / £ e (\x\,t)dp(t,x) = 1, which concludes the 
lemma by letting e go to 0. 

□ 

From this lemma, we have: 

/ Vy?(t,x) -x 1 - dp n (dt,x) + / Vip(t,V$[p n (t)](x)) ■ or 1 dtdx 

i[0,T]xM 3 J[0,T]xQ. 

< C(T)\\ip\\ L iQp jT \ } ci(B R(X) y 

Thus from Definition 12.11 equation © we know that for any time t > 0, 
d t p n (t, •) is bounded in the dual of -^ 1 ([0, T], C 1 (1R 3 )) and thus in the dual of 
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L 1 ([0,T], W 2,P (M. 3 )) for p > 3 by Sobolev embeddings. Thus for some p* > 1 
we have 

9 t p n eL°°([0,T],W- 2 ' p '(R 2 )). 

With the two above results, and using classical arguments of functional 
analysis (see we can obtain the following lemma: 

Lemma 2.4. Let the sequence (p n ) ng N be as above, there exists p G C([0, T],V— 
w*) and a subsequence (p nk )keN, such that for allt G [0,T], p nk (t) converges 
to p{t) in the weak-* topology of measures. 

With this lemma, we need to show that for all <p G C£°([0, T] x M 3 ) 
we have V<p(t, V$[p n (t)]) converging to V<p(t, V&[p(t)]) whenever p n (t) con- 
verges weakly-* to p(t). This last step will be a consequence of the following 
stability theorem: 

Theorem 2.5 (Brenier, |3J). Let Q be as above. Let (p n ) n£ N be a sequence 
of probability measures on M. d , such that Vn, J(l + \x\ 2 )dp n < C, let $ n = 
$[p n ] and \l/ n = \I/[p„] be as in Theorem li.il If for any f G C ,0 (lR d ) such 
that \f(x)\ < C(l + |x| 2 ), f fp n — > f pf, then $„ — > $[p] uniformly on each 
compact set of Q and strongly in M. d ), and — * ^[p] uniformly on 

each compact set ofM. d and strongly in W^' c (M d ). 

From this result, we obtain that the sequence V^fpn] converges strongly 
in and almost everywhere (because of the convexity of $[/?]) to V$[p]. 

Thus V<p(£, V$[p n ]) converges to Vip(t, V$[p]) in L 1 (f2) and one can pass to 
the limit in the formulation of Definition 12.11 This ends the proof of point 2 
of Theorem O 

Existence of solutions 

We show briefly the existence of a solution to the Cauchy problem in the 
sense of Definition 12.11 Indeed given p° the initial data for the problem that 
we want to solve, by smoothing p°, we can take a sequence p° of initial 
data belonging to L^R 2 ), uniformly compactly supported and converging 
weakly-* to p°. We know already from [T], [TU], [HI that for every p°, one 
can build a global weak solution of (0111 EJ), that will be uniformly compactly 
supported on [0,T] for all T > 0. This sequence will also be solution in the 
sense of Definition ^. 11 We then use the stability result, and conclude that, up 
to extraction of a subsequence, the sequence p n converges in C([0, T], V — w*) 
to a weak measure solution of SG with initial data p°. This achieves the proof 

of Theorem O □ 
Remark: One can prove in fact the more general result, valid for non 
linear functionals: 
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Proposition 2.6. Let F e C°{QxR d ), such that \F(x,y)\ < C{l + \y\ 2 ), let 
(Pn)neN be a bounded sequence of probability measures, Lebesgue integrable, 
with finite second moment. Let p be a probability measure with finite sec- 
ond moment, such that for all f G C°(M ) such that \f(x)\ < C(l + \x\ 2 ), 



3 Continuous solutions 

What initial regularity is necessary in order to guarantee that the velocity 
fields remains Lipschitz, or that the flow remains continuous, at least for a 
short time ? The celebrated Youdovich's Theorem for the Euler incompress- 
ible equation shows that when d = 2, if the initial vorticity data is bounded 
in L°°, the flow is Holder continuous, with Holder index decreasing to as 
time goes to infinity. This proof relies on the following regularity property 
of the Poisson equation: if A0 is bounded in L°°, then V0 is Log-Lipschitz. 
This continuity is enough to define a Holder continuous flow for the vector 
field V0 -1 . Such a result is not valid for the Monge- Ampere equation. As far 
as we know, the optimal regularity result for Monge- Ampere equations is the 
following: 

3.1 Regularity of solutions to Monge- Ampere equation 
with Dini-continuous right hand side 

Theorem 3.1 (Wang, |22j). Let u be a strictly convex Alexandrov solution 



with p strictly positive. Ifw(r), the modulus of continuity of p, satisfies 





of 



det D 2 u = p 



(7) 




8 



then u is in C? 
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We will work here in the periodic case. In this case, u the solution of 
(J2j) will be ty[p] of Theorem 11.21 The arguments of [7j, [H], adapted to the 
periodic case, show that ty[p\ is indeed a strictly convex Alexandrov solution 
of solution of Therefore we obtain the following corollary of Theorem 

ECU 

Corollary 3.2. Let p e V(T d ) be such that 

< m < p < M, 

f 1 w ( r ) j ^ 

/ —^-dr = C < oo. 

Jo r 

where m,M,C are positive constants. Let be as in Theorem \l. 6 A We 
have, for some constant Ti depending on m, M, C 

\Mp]\\c*(Td) < n. 



3.2 Result 

We will now prove the following: 

Theorem 3.3. Let p° be a probability on T 3 , such that p is strictly positive 
and satisfies the continuity condition (0). Then there exists T > and Ci, C2 
depending on p° , such that on [0, T] there exists a solution p(t,x) of SG that 
satisfies for all t G [0, T] : 

J^v^rldrKCt, \m,-)\\cv>)<c 2 , 

where w(t,r) is the modulus of continuity (in space) of p(t,.). 



Proof of Theorem 13.31 

Let us first sketch the proof: If \I/ <G C 2 , then the flow t — > X(t, x) generated 
by the velocity field [V^a;) — x} 1 - is Lipschitz in space. Since the flow is 
incompressible, we have p(t,x) = p°(X -1 (t, x)). 

Now we use the following property: If two functions /, g have modulus of 
continuity respectively Wf,w g then g o f has modulus w g o wj. 
Thus if is Lipschitz, we have w p o x~i(t) ^ w p o(L •) with L the Lipschitz 

constant of X _1 (t) and condition (JHJ) remains satisfied. 

Remark 1: Note that Holder continuous functions satisfy the condition 

©• 

Remark 2: Note also that we do not need any integrability on Vp and the 
solution of the Eulerian system still has to be understood in the distributional 
sense. 
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A fixed point argument 



Let us introduce the semi-norm 



i 




/"lie 



dr 



(9) 



r 



defined on T(T 3 ), where we recall that w p is the modulus of continuity of p. 
We denote Vc the set V equipped with this semi-norm, i.e. 



From now, we fix p° a probability density in Vc, satisfying m < p° < M, 
where m and M are strictly positive constants. Let p be a time dependent 
probability density in L°°([0, T]; Vc), such that m < p(t) < M for all t, we 
consider the solution p of the initial value problem: 



From Theorem 13 . 1 1 and its corollary, the vector field v|/t] = (V^ [p](x) — x) ± is 
C 1 uniformly in time, therefore there exists a unique solution to this equation, 
by Cauchy-Lipschitz Theorem. This solution can be built by the method of 
characteristics as follows: Consider the flow X(t,x) of the vector field v[/x], 
then p{t) is p° pushed forward by X(t), i.e. pit) = p° o X _1 (t). From the 
incompressibility of v[/x] the condition m < p° < M implies that for all 
t E [0,T], m < p(t) < M. 

The initial data p° being fixed, the map p i— > p will be denoted by T. 

The spatial derivative of X, D X X satisfies 



and the same bound holds for X(t) 1 . Since Wf og < Wf o w g , and writing 
C t = exp(tsup sg[ot] \D x v{p]\), we obtain w p{t) (-) < w p o(C r ), and 



V c = {p E P(T 3 ), ||p|| c < oo}. 



d t p+(V^[p]{x)~x) ± -Vp = 0, 
p(t = 0) = P °. 



(10) 
(11) 



d t D x X = D x v[p}(X)D x X, 



therefore we have 



D x X(t)\ < exp(t sup \D x v[p]{s)\) 



s<E[0,t] 




1). 
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(using that Vr, w p (r) < M — m). Therefore, 

\\p(t)\\c<\\p°\\c + (M-m)(C t -l). 

Now from Corollary 13. 21 and m, M being fixed, there exists a non-decreasing 
function H such that 

H/i]|| C i <H(\\fi\\e), 

and so C t < exp(tTi,(\\ii\\L°°([o,t];Vc)))- Hence we can chose Q > 1, and then 
T such that 

|| p 0|| c + (M _ m) ( exp(T H(Q\\P°\\c)) " 1) = QIIpIc 

Note that for Q > 1, we necessarily have T > 0. Then the map JF : p i— > p 
goes now from 

.A = {//, ||/i|U-([o,T];P c) < Q||p°||e, m < fi < M} 

into 

B = {p, ||p(*)||c < ||p°||c + (M - m) (exp(t W(Q||p°|| c )) - 1 ) , Vt G [0, T] } , 

and with our choice of T = T(Q), we have B C A. Moreover from the 
unconditional bounds 

p<M, 

llvNII L°°([0,T]xT 3 ) 

(see the remark after Theorem II .21 for the second bound) and using equation 
f]"0|) . we have also ||0tp||i,°°([o,T];W<'- 1 ' <x> ) < K(M) whenever p = J-'(p). 

Call ^4 (resp. B) the set „4 fl {p, \\d t p\\ l°° ([o^w- 1 ' 00 ) < ^X^0}> (resp. 
B fl {p, HQt/oHioo^o^^-i.oo) < K{M)})\ we claim that 

• F(A) dBdA, 

• A is convex and compact for the C°([0, T] x T 3 ) topology, 

• T is continuous for this topology, 

so that we can apply the Schauder fixed point Theorem. We only check the 
last point, the second being a classical result of functional analysis. So let 
us consider a sequence (p, n )neN converging to p G A, and the corresponding 
sequence (p n = J r (p n ))neN- The sequence p n is pre-compact in C°([0, T] x T 3 ), 
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from the previous point, and we see (with the stability Theorem 12. 5|) that it 
converges to a solution p of 

%> + V-(pv[/i]) = 0. 

But, v[/i] being Lipschitz, this solution is unique, and therefore J 7 (/J ln ) con- 
verges to J-'(fi), which proves the continuity of J 7 , and ends the proof of 
existence by the Schauder fixed point Theorem. 

□ 

We state here without proof some consequences of the previous result: 
Corollary 3.4. Let p° G V(T 3 ), such that < m < p < M. 

1. If p° G C a ,a G]0, 1], for T* > depending on p° , a solution p(t,x) to 
J3jp) exists in L°°([0,T*[,C Q (T 3 )). 

2. If p° G W 1,p ,p > 3, for T* > depending on p°, a solution p(t,x) to 

exists in L°°([0,T[, W 1,P (T 3 )). 

3. If p° e C k ' a ,a G]0,1],A; G N, forT* > depending on p° , a solution 
p(t,x) to exists in L°°([0, T% C k > a (T 3 )). 

Moreover, for these solutions, the velocity field is respectively in C 1,a (T 3 ), 
W 2 ' P (T 3 ), andC k+1 ' a (T 3 ) on [0,T*[. 

4 Uniqueness of solutions to SG with Holder 
continuous densities 

4.1 Result 

Here we prove the following theorem: 

Theorem 4.1. Suppose that p° G ^(T 3 ) with < m < p° < M, and belongs 
to C a (T 3 ) for some a > 0. From Theorem \3.3i for some T > there exists 
a solution p to SG in L°°([0, T], C a (T 3 )). Then every solution of SG in 
L°°([0, T'], ^(T 3 )) forT' > 0,/3 > with same initial data coincides with p 
on [0, inf {T,T'}]. 

Remark 1: The uniqueness of weak solutions is still an open question. 

Remark 2: Our proof of uniqueness is thus valid in a smaller class of 
solutions than the one found in the previous section, the reason is the follow- 
ing: during the course of the proof, we will need to solve a Monge- Ampere 
equation, whose right-hand side is a function of the second derivatives of the 
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solution of another Monge- Ampere equation. In Theorem 13 .11 if u is solution 
to (JJJ) with a right hand side satisfying (jHJ), although u G C 2 , it is not clear 
that the second derivatives of u satisfy (jHJ)- Actually, it is even known to be 
wrong in the case of the Laplacian (for a precise discussion on the subject, 
the reader may refer to (H|). However, from Theorem 14.31 below, if p G C a 
then u G C 2 ' a . 

What we actually need is a continuity condition on the right hand side of 
(JTj) such that the second derivative of the solution u satisfies (jHJ)- This may 
be a weaker condition than Holder continuity, however the proof would not 
be affected, therefore it is enough to give it under the present form. 

Proof of Theorem 14.11 

Let pi and p 2 be two solutions of (0 H EJ) , in L°°([0, T], C^T 3 )) that coincide 
at time 0. Let Xi,X 2 be the two corresponding Lagrangian solutions, (i.e. 
solutions of (ffP ). The velocity field being C 1 , for all t G [0,T], X^t, •) and 
X 2 (t, •) are both C 1 diffeomorphisms of T d . 

We call v x (resp. v 2 ) the velocity field associated to Xi (resp. X 2 ), 
Vi(t, x) = x) - x} ± ,i = 1, 2. We have 



We want to obtain a Gronwall type inequality for ||Xi — X 2 ||i2. Since vi 
is uniformly Lipschitz in space (from Theorem 13. 3|) . the first bracket is esti- 
mated in L 2 norm by C||Xi — X 2 || L 2. 

We now need to estimate the second term. We first have that 



and since p 2 is bounded, we need to estimate ||V^i — V^IU 2 - This will be 
done in the following Proposition: 

Proposition 4.2. Let Xi,X 2 be mappings from T d into itself, such that the 
densities pi = X^dx.i = 1,2 are in C a (Y d ) for some a > 0, and satisfy 
< m < pi < M. Let \l>i, i = 1, 2 be convex such that 



dt(Xi — X 2 ) 



vi(A 1 )-v 2 (A 2 ) 

(v^) - Vl (X 2 )) + ( Vl (A 2 ) - v 2 (X 2 )). 




detD 2 ^ = pi 
in the sense of Theorem M . 11 i.e. = ^[pj]. Then 



V* 1 -V*2\\i?<C\\X 1 -X 2 \\v, 



where C depends on a (the Holder index of pi), ||pi||c<*(T d )j m an d M. 
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Before giving a proof of this result, we conclude the proof of the Theorem 
14.11 The Proposition 14.21 implies immediately that 

\\a t (x 1 -x 2 )\\v<c\\x 1 -x 2 \\v, 

and we conclude the proof of the Theorem by a standard Gronwall lemma. 

□ 

4.2 Energy estimates along Wasserstein geodesies: Proof 
of Proposition 14. 2L 

In the proof of this result we will need the following result on optimal trans- 
portation of measures by gradient of convex functions: 

Theorem 4.3 (Brenier, |3J, McCann, |2()| . Cordero-Erausquin, |9j, 
Caffarelli,|6j). Let p 1; p 2 be two probability measures on T d , such that p\ is 
absolutely continuous with respect to the Lebesgue measure. 

1. There exists a unique up to a constant convex function <p such that 
(f) — | ■ | 2 /2 is Z d periodic, satisfying V0#pi = p 2 . 

2. The map V0 is the solution of the minimization problem 

inf / pi(x)\T(x) — x\^ d dx, (12) 

T#pi=p 2 Jjd 

and for all x G M d ; | V0(x) — x\ T d = | V0(x) — x\ R d. 

3. If p±, p 2 are strictly positive and belong to C a (T d ) for some a > then 
<fi G C 2 ' a (Y d ) and satisfies pointwise 

p 2 (V0) det£> 2 = pi. 

For complete references on the optimal transportation problem (|T2| and 
its applications, the reader can refer to [2Tj . 

Remark 1: the expression \- \ jd denotes the Riemannian distance on the 
flat torus, whereas | • \ R d is the Euclidian distance on W d . The second assertion 
of point 2 means that, for all x G M d , |V0(x) — x\ < diam(T d ) = \fd/2. 

Remark 2: Here again, note that since — | ■ | 2 /2 is periodic, the map 
x I— > V0(x) is compatible with the equivalence classes of M d /Z d , and therefore 
is defined without ambiguity on T d . 
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Wasserstein geodesies between probability measures 

In this part we use results from [2], PH]. Using Theorem 14.31 we consider the 
unique (up to a constant) convex potential <fi such that 

- | • | 2 /2 is Z d - periodic. 
We consider, for 9 G [1, 2], <f)g defined by 

I 1 2 

^ = (2-0)^- + (0-l)«/>. 

We also consider, for 6* G [1, 2], p 9 defined by 

p e = V0 e #pi. 

Then p# interpolates between pi and p 2 . This interpolation has been in- 
troduced in [2| and JH] as the time continuous formulation of the Monge- 
Kantorovitch mass transfer. In this construction, a velocity field v g is defined 
dp e a.e. as follows: 

V/ G C°(T d ; R d ), J Pe v e ■ f = J Pl f(V<p e ) • dgVfo- (13) 

It is easily checked that the pair pg,v$ satisfies 

d e pe + V • (pev e ) = 0, 
and for any 9 G [1, 2], we have (see [2]): 

^/ Pe\ve\ 2 = lf Pi|V0(x) -x| 2 = Wf(pi,p 2 ), 

Z J^d Z J-jd 

where ^2(^1,^2) is the Wasserstein distance between pi and p 2 , defined by 
W*{ Pl ,p 2 ) = inf ( / pi(x)|T(x) -x| 2 



The Wasserstein distance can also be formulated as follows: 



Wi( PuP2 )= mi^j^-Y^A 



where the infimum is performed over all maps Y U Y 2 : T d i-> T d such that 
Yi#dx = pi, i = 1,2. From this definition we have easily 



W^(pi,p 2 ) < J \X 2 (t,a)-X l (t,a)\ 2 da, 
and it follows that, for every 9 G [1,2], 



/ p e \ve\ 2 = W 2 2 {p 1 ,p 2 ) < \\X 2 -X 1 \\ L 2. (14) 
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Regularity of the interpolant measure pg 

From Theorem 14.31 for pi, p 2 G C@ and pinched between the positive postive 
constants m and M , we know that <t> G C 2,13 and satisfies 



det D 2 c 



Pi 



P2W) 



We now estimate pg = pi[det D 2 cj)g] 1 . From the concavity of log(det(-)) 
on symmetric positive matrices, we have 

det D 2 (j)g = det((2 - 9)1 + (9 - l)D 2 <f)) 

> [det D 2 ^} 9 - 1 
m 

> — . 

- M 

Moreover, since <fi G C 2 , det D 2 cf)g is bounded by above. Thus pg is uniformly 
bounded away from and infinity, and uniformly Holder continuous. 

Final energy estimate 

If we consider, for every 9 G [1, 2], ^>g solution of 

detD 2 ^g = p e , (15) 

in the sense of Theorem II. 2\ then tyg interpolates between and \&2, and 
tyg G C 2 '^ uniformly, from the regularity of pg. We will estimate dgV^g by 
differentiating (fl"5|) with respect to 9: for M, iV two d x d matrices, t G M, 
we recall that 

det(M + tN) = det M + t (trace M l co N) + o(t), 

where M co is the co-matrix (or matrix of cofactors) of M. Moreover, for any 
/ G C 2 (IR' i ; R), if M is the co-matrix of D 2 f, it is a common fact that 

d 

Vj e[l..d],J2diMij = 0. (16) 



i=l 



Hence, denoting Mg the co-matrix of $j, we obtain that dg^/g satisfies 

V-(MgVdg^g) = dg P g(t) 

= -V-(pgVg). (17) 
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From the C 2 ^ regularity of tyg, D 2 ^g is a C 13 smooth, positive definite matrix, 
and its co-matrix as well. Thus the problem (JT7|) is uniformly elliptic. If we 
multiply by dg^ffg, and integrate by parts we obtain 

V%V e M e Vde^e = - Jvd e ^ e - v ePe . 

Using that Mg > XI for some A > 0, and combining with the inequality (fT4*j) 
above, we obtain 

\\Vdg^g(t)\\ L 2 < A _1 ||p^|| L 2 

1/2 



< A 1 \\X 2 - Xi\\ L 2 fsup UpelliooJ 



The constant A -1 depends on m, M, /3, {||pi||c' 3 > i = 1> 2}, and is thus bounded 
under our present assumptions. We have already seen that pg is uniformly 
bounded, and we finally obtain that 

||Wi-W 2 ||l2 <C\\X 1 -X 2 \\ L 2, (18) 

this ends the proof of Proposition 14.21 



□ . 

Remark 1. In ^7], the author obtains also (weaker) estimates of the type 
of Proposition 14. 2| for discontinuous densities pi,p 2 . 

5 Uniqueness of solutions to the 2-d Euler 
equations with bounded vorticity: a new 
proof 

This proof adapts easily to the case of 2-d Euler equation with bounded 
vorticity, giving a new proof of the uniqueness part in Youdovich's theorem. 
We start now from the following system: 

d t p + W ± • Vp = 0, (19) 
P = AV, (20) 
pit = 0) = p°. (21) 

For simplicity, we restrict ourselves to the periodic case, i.e. x G T d , p,ip 
periodic, this implies that p has total mass equal to 0. We reprove the 
following classical result: 

Theorem 5.1 (Youdovich, |23j). Given an initial data p° £ L°°(T 2 ) sat- 
isfying J T2 p° = 0, there exists a unique solution to ( Tlffi W(\ H7|) such that p 
belongs to L% C (R + x T 2 ). 
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Proof of Theorem 15.11 

We consider two solutions pi, ipi and p 2 ,ip 2 , such that p iy i = 1,2 are bounded 
in L°°([0,T] x T d ). In this case the velocity fields v$ = both satisfy 

V(x,?/)gT 2 , \x-y\<\, \vi(x) - Vi(y)\ < C\x-y\\og- — - — -. 

2 \x — y\ 

This implies that the flows (t,x) i— > Xi(t,x) associated to the velocity fields 
Vi = Vip^ are Holder continuous, and measure preserving. Moreover, one 
has, for all t G [0,T], Pi {t) = Xi(t)#p°. 

Applying the same technique as before, we need to estimate ||VV>i — 
V^2||l 2 (t 2 ) in terms of ||Xi — -X"2||x, 2 fx 2 )> 111 the present case, the energy 
estimate of Proposition 14.21 will hold under the weaker assumptions that the 
two densities are bounded. 

Proposition 5.2. Let Xi,X2 be mappings from T d into itself, let p° be a 
bounded measure with a density in L°° with respect to the Lebsgue measure, 
and with f Td p° = 0. Let pi = X i# p°,i = 1,2. Let ip^i = 1,2 be periodic 
solutions of Aip i = pi,i = 1,2, then we have 

|| W>i - VV>2||z2(t<*) < (2||p°|U°c max-nipxllz^, ||p 2 ||l-}) 1/2 ||-Xi - X 2 \\ L 2 (J d ) . 

Remark: In other words, this proposition shows that for p\,p2 bounded, 
the if -1 norm of p\ — p 2 is controlled by some 'generalized' (since here we 
have unsigned measures) Wasserstein distance between pi and p 2 . 

To conclude the proof of Theorem 15. 1[ note first that for all C > 0, we 
can take T small enough so that ||X 2 — Xi||£,oo([ 0i t]xt 2 ) < C. Now we have 
for the difference X\ — X 2 , as long as \X\ — X 2 \ < 1/2, 

\\d t (X 1 -X 2 )\\ L 2 

< \\vm x i)-vM x *)\\v + liv^i(x 2 ) - v^ 2 (x 2 )|| L2 

< d|| \X, - X 2 \ logdXi - x 2 \)\\ L , + c.wx, - X 2 \\ L t, 

where, to evaluate the second term of the second line, we have used the fact 
that 

||vVi(x 2 )-vv 2 (x 2 )|U 2 = ||Wi - 

and then applied Proposition 15.21 

We just need to evaluate || \X% — X 2 \ log(|Xi — X 2 |)||l2. We take T small 
enough so that ||X 2 — Xi|| iO o([ ,T]xT 2 ) < l/ e and notice that x 1— > xlog 2 x is 
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concave for < x < 1/e, therefore by Jensen's inequality we have 

f \X 2 -X 2 \ 2 \og\\X 1 -X 2 \) 

and some elementary computations finally yield 

d t \\X 2 -X 1 \\ L 2<C\\X 2 -X 1 \\ L . log ||y l —. . 

11^2 - -m||l 2 

The conclusion Xi = X 2 follows then by standard arguments. 

5.1 Energy estimates along Wasserstein geodesic: Proof 
of Proposition 15.21 

The proof of this proposition is very close to the proof of Proposition 14.21 
and we will only sketch it, insisting on the specific points. Here the densities 
Pi can not be of constant sign, since their mean value is zero, hence we 
introduce p°' + (resp. p 0, ~) the positive (resp. negative) part of p°. Then we 
introduce pf = Xj#p 0,± . Note that if the mappings Xj were injective, (which 
is the case in our present situation) we would have pf that coincides with 
the positive/negative parts of p,, but this can be wrong if Xj is not injective. 
However what remains is that pi = pf — pi . Now, pf, i = 1,2 are 4 positive 
measures of total mass equal to say M, with M < oo. 

Wasserstein geodesic 

We interpolate between the positive parts pf , and the negative part is han- 
dled in the same way. As before we introduce the density pf{t) that inter- 
polates between pf(t) and pf(t). In this interpolation, we consider vf such 
that 

depf + V • (pfvf) = 0, 

and we introduce as well pg,Vg. Then pe = pf — Pg has mean value 0. Let 
the potential ifjQ be solution to 

AVe = p e - (22) 

Note that pe has mean value zero therefore this equation is well posed on T 2 , 
moreover ipg interpolates between tpi and ip 2 - 
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Bound on the interpolant measure p# 

Instead of interpolating between two smooth densities, we interpolate be- 
tween bounded densities, and use the following result from [TH] : 



Proposition 5.3 (McCann, [19]). Let pj~ be the Wasserstein geodesic link- 
ing Pi to p2 defined above. Then, for all 9 G [1,2], 

\\p~e < max{||p^|| L ^, 11/4 IU°°} • 
The same holds for p~[ , p$ . 

Remark: This property is often referred to as 'displacement convexity'. 

Energy estimates 

Now by differentiating (}22|) with respect to 6, we obtain 

Adetpg = d e pe = -V • (pjv£ - p e v e ), (23) 



with the interpolating velocity defined as in (fT3|) . and satisfying for all 

Be [1,2], 



Multiplying (J23j) by dgipe, and integrating over 9 G [1,2], we obtain 

|| V^l -V^ 2 1| L2 (T d) < / ||PgV|| £ 2 + \\PgVg\\ L 2 

J 0=1 



< W 2 (pf,pt) (suv\\pj\\ L o?J 

+ W 2 (pi,p 2 ) ^sup Hp^Hl-^ • 

Note that the energy estimate is easier here than in the Monge- Ampere case, 
since the problem is immediately uniformly elliptic. 

The mappings Xi satisfy Xj#p = p^ and Xj#(pg ) = pf. Hence, 

wM.pt) < Jpt\Xi-x 2 f. 

Using Proposition 15. 3[ we conclude: 
IIWi-V^H L 2 (T d ) 
< 2||X 2 -Xi|| L 2^||p || L o O max{||pi|| L oc, ||p 2 |U°°} 
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This ends the proof of Proposition 15.21 Note that in our specific case, Xi 
are Lebesgue measure preserving invertible mappings, therefore \\pf = 
\\Po\\l°°, and the estimate can be simplified in 



Remark: This technique can be used to conclude uniqueness for many 
non-linear systems, where a transport equation and an elliptic equation are 
coupled. The velocity field is the gradient of a potential satisfying a elliptic 
equation whose right hand side depends smoothly on the density. For exam- 
ple, we have uniqueness of solutions to the Vlasov-Poisson system and Euler- 
Poisson system with bounded density and bounded velocity. The Vlasov- 
Monge- Ampere and Euler-Monge- Ampere systems have also been studied by 
the author (J3], and the same technique apply to yield uniqueness for 

solutions with C a density and bounded velocity. Note however that to en- 
force uniform ellipticity, we need for the Monge- Ampere equation the density 
to be bounded by below which is not the case for the Poisson equation. 

6 Convergence to the Euler equation 
6.1 Scaling of the system 

Here we present a rescaled version of the 2-d SG system and some formal 
arguments to motivate the next convergence results. Here x G T 2 ,t G 1R + 
and for v = (vi,v 2 ) G M 2 , v 1 - now means (— V2,V\). Introducing ip[p] = 
*&[p] — \x\ 2 /2, where ty[p] is given by Theorem 11.21 the periodic 2-d SG 
system now reads 



If p is close to one then if) should be small, and therefore one may consider 
the linearization det (I + D 2 ip) = 1 + Aip + 0(\D 2 ip\ 2 ), that yields Aip ~ p—1. 
Thus for small initial data, i.e. p° — 1 small, one expects if), p = p — 1 to stay 
close to a solution of the Euler incompressible equation EI 



\\Vlf>l — V^2||L2( T d) < 2||po||l°° II X 2 — Xi\\ L 2(jdy 



□ 



d t p + V • (pV^ 1 ) = 0, 
det(J + DV) = P- 



% + V-(pV0 ± ) = O, 
A0 = p. 



(24) 
(25) 
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We shall rescale the equation, in order to consider quantities of order one. 
We introduce the new unknown 

p e (t,x) = -(p(-,x) - 1), 
e e 

^j e (t,x) = -if){-,x). 



Then we have 



and we define 4> e by 



so that 



p(t) = l + ep e (et), 
y[p](t) = \x\ 2 /2 + er(et), 

e0 e = |x| 2 /2-$[p], 



V0 £ = V^(V<%]). (26) 

Hence, at a point x G T 2 , V0 €_L is the velocity of the associated dual point 
V$[p](o;). The evolution of this quantities is then governed by the system 
SG € 

d t p e + V ■ (p e W> £± ) = 0, (27) 
det(/ + eDV) = 1 + ep e . (28) 

Remark: Note that this system admits global weak solutions with initial 
data any bounded measure p e0 , as long as 

P e ° = 0, (29) 
P eu > --• (30) 



Note also that if the pair (p, (f>) is solution to the EI system fl2"^ l23j) . so 

1 ,t , l„t 



is the pair ( -p(-,x), —(/)(-, x 
\e e e e 

We now present the convergence results. We show that solutions of SG e 
converge to solutions of EI in the following sense: if p e0 , the initial data of 
SG t , is close (in some sense depending on the type of convergence we wish 
to show) to a smooth initial data p° for EI, then p e and p remain close for 
some time. This time goes to oo when e goes to 0. 

We present two different versions of this result: the first one is for weak 
solutions of SG e , and the second one is for Lipschitz solutions. 
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6.2 Convergence of weak solutions 

Theorem 6.1. Let (p t ,ip t ) be a weak solution of the SG e system Wty) . 
Let (p,0) be a smooth C 3 ([0, T] x T 2 ) solution of the EI system |§H|. Let 
(p e be obtained from as in [2b)) . let H e (t) be defined by 

He{t) = \ [ |V0 e -V0| 2 , 
- JT 2 

then 

H e {t) < {H t (0) + Ce 2/3 (l+t))expCt 
where C depends on sup 0<s<t {\\D 3 (j)(s), D 2 d t (p(s) ||l°o(t 2 )}- 

Remark 1: Note that V0 e-L (t, x) is the velocity at point V$[p] = x — eV0 e . 
Thus we compare the SG e velocity at point x— eV0 e (the dual point of x) with 
the EI velocity at point x. Our result allows also to compare the velocities 
at the same point, by noticing that 

G e (t) = 



< 

using the smoothness of 4>, and if v S9e , v ei are the respective velocities of the 

SG e and EI systems, G e = p e \^ sg , — v ei | 2 . 

Remark 2: The expansion det(7 + D 2 vp) = 1 + + 0(\D 2 i>\ 2 ), used 
above to justify the convergence relies a priori on the control of D 2 ijj in the 
sup norm. But in the Theorem I6.1| the initial data must satisfy Vip e close 
in L 2 norm to a smooth divergence free velocity: this condition means that 
D 2 ip e is close in H _1 norm to D 2 <ft, which is smooth. This control does not 
allow to justify the expansion det(J + D 2 ip) = 1 + + 0(\D 2 ip\ 2 ), but we 
see that the result remains valid. 

Proof of Theorem 16.11 

In all the proof, we use C to denote any quantity that depends only on (f>. 
We use the conservation of the energy of the SG t system, given by 

E(t) = [ |V0f . 
Jt 2 



\ / P |V^-V0| 2 
^ Jt 2 

1 [ |V0 e - V0(x-eV0 e ) 

2 Jt 2 

C(H e (t) + e 2 ) 
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This fact, although formally easily justified, is actually not so straightforward 
for weak solutions, and has been proved by F. Otto in an unpublished work. 
The argument is explained in p3|. Therefore E(t) = Eq. The energy of the 
smooth solution of EI is given by 



|V# 



2 



and also conserved. For all smooth 9, we will use the notation: 



< D 2 6 > (t, x) = I (1 - s)D 2 9(t, x - seV0 e (t, x)). 

's=0 



Thus we have the identity 



p'9 = / 6(x-eV(j) e ) (31) 

T 2 JT 2 

= I 9 - e [ V0 ■ V0 e + e 2 I < D 2 9 > V0 e V0 e . (32) 

Using the energy bound, the last term is bounded by e 2 \\D 2 9\\ L o C ^E . Then 
we write 



Using the identity (j32j), we have for all smooth 9, 



e W ■ V0 e = - / p e 9+ 9 + e z < D z 9 > 

./' ' 2 JT 2 JT 2 J •■' 

hence, replacing 9 by in this identity, we get 

±H € {t) = ~f [p £ 0-0-e 2 < J D 2 0>V0W]- 
dt e dt Jj2 

We can suppose without loss of generality that / <j){t,x) dx = 0. Then if 

Jt 2 

we define 

Q e {t) = [ e<D 2 0> V0 £ V0 e , 
Jt 2 

(note that |Q e (*)| < Ce), we have 
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Hence we are left to compute 

e at j T 2 e J T 2 

1 



p e W e± ■ V0 - eV0 e • Vd t cj) + e 2 < D d t <f) > V0V0 

e ./t 2 



1 



p e V^ eL -V(j)- / V(j) e -Vd t (f) + 0(e) 
e Jt 2 jt 2 

= T 1 + T 2 + 0{e), 

where at the second line we have used (|27|) for the first term and (|3*2*j) with 
6 = d t (f) for the second and third term. (Remember also that we assume 
fdt$ = 0.) 

We will now use the other formulation of the Euler equation: v = 
satisfies 

dtV + v • Vv = — Vp. 
After a rotation of tt/2, this equation becomes: 

d t V0 + D 2 0V0 L = V P ± , 

thus for T2 we have 

Jt 2 



v</- en2 

'T 2 

For Ti, using (|2*^jl and (J52)l . we have 



/t 2 



/ W> e± (x-eV0 e ) • V0(x-eV0 e ) 
Jt 2 



T 2 



where A is defined by 



A = J Vcf> e± (^D 2 (j) - jf J D 2 0(s - seV0 e ) cfe ^ V0 e . (33) 

The term J T2 V0 e_L ■ V0 vanishes identically. Concerning A, we claim the 
following estimate: 
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Lemma 6.2. Let A be defined by Uty) . then 

|A| < C(e§ + H e ), 
where C depends on ||-D 3 <^||loo . 

We postpone the proof of this lemma after the proof of Theorem 16.11 We 
now obtain 

^-(H e (t) + Q e {t)) < [ (V^ - V<p e± )D 2 4>V<p e + CH e + Ce 2 '\ 
at J T 2 

Noticing that for every 9 : T 2 i— > R we have 

/ V0 ± D 2 p7j> = [ V6 ± - V(^|V0| 2 ) = 0, 
Jt 2 Jt 2 2 

we find that 

/ (V0 X - V<p e± )D 2 0\70 e = [ (V0 X - V0 e± )Z^ 2 0(V0 e - V0), 
Jt 2 Jt 2 

hence 

^-(H e (t) + Q £ (t)) < -f (V0 e± - V0 L )D 2 0(V(j) e - V0) + + Ce 2 ' 3 
at j T 2 

< C (H e (t) + Q e (t) + e 2 / 3 ) 
using that < Ce. Therefore 

H e {t) + Q £ (t) < (F e (0) + Q £ (0) + Ce 2/3 t) exp(Ct) 

and finally 

# £ (t) < (H £ (0) + Ce 2 / 3 (1 + t)) exp(Ct) 



and the result follows. Check that the constant C depends only on 
sup 0<s<t {\\D 3 (j), D 2 d t (j)\\L°°(T 2 )}- This ends the proof of Theorem 16.1 



□ 



Proof of Lemma 16.21 

First we show that if Q(K) = / |V0 e | 2 , then 

Q(R) <C [ |V0 e - V0| 2 + ^. (34) 
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Indeed, / |V0 e | 2 < C, implies that meas{|V> e | > R} < C— . Since \V(p{t,x)\ < 
C for (t, x) e [0, V] x T d , we have 

Q{R) < I |V0| 2 + / |V0 e -V0| 2 

-/{|Vai e |>.R} ^{|V</> E |>,R} 

y iv^-v0i 2 . 



c 

< — 
- R 2 



Hence plj) is proved. 
Then, letting 



K{x) = D 2 <j) - [ D 2 <f)(x - seV<f) € ) ds, 

Js=0 



we have 



A<CB(R)+ [ \K{x)\\V<f) e \ 2 

J\W^\<R 

with \K{x)\ < Ce\V(j) € \ thus 

A < Ce [ \V4> e \ 3 + CQ(R) 

J\V^\<R 

<c(eRj |V0f + + J W - V0| 2 
<c(eR + ^ + J |V0 £ -V0| 2 ) 
for all R, so for R = e -1 / 3 we obtain: 

A < Ce 2/3 + C J |V0 e - V0| 2 . 



This proves Lemma |6l 

□ 



6.3 Convergence of strong solutions 

We present here another proof of convergence, that holds for stronger norms. 
Let us consider as above the solution (p, cj)) to Euler: 

d t p + V-(pV^) = 0, 
A0 = p, 
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and we recall the SG £ system 



d t p e + V ■ (p e W> £± ) = 0, 
det(I + e£>V) = 1 + ep e - 



We have then 



Theorem 6.3. Let (p, 0) be a solution of EI, such that that p £ Cf oc (R + x 
T 2 ). Let p e0 be a sequence of initial data for SG e satisfying Wfy and 
p e0 -p° 

such that is bounded in W l, °°{T 2 ). Then there exists a sequence 

€ 

(p e ,ip e ) of solutions to SG e that satisfies: for all T > 0, there exists €t > 0, 
such that the sequence 

e e 

for < e < tT is uniformly bounded in L°°([0,T], W /1, °°(T 2 )). 

Remark: In the previous theorem, we obtained estimates in L 2 norm, here 
we obtain estimates in Lipschitz norm. Estimates of higher derivatives follow 
in the same way. 

Proof of Theorem 16.31 

We expand the solution of SG £ as the solution of EI plus a small perturbation 
of order e and show that this perturbation remains bounded in large norms 
(at least Lipschitz). We first remark the the assumption on p implies that 
VT > 0, G L°°([0,T];C 3 (T 2 )). Let us write 

P e = P + epi 
ip e = 4> + eipi- 

Rewritten in terms of Pi,ipi, the SG t system reads: 

d t pi + (V0 + eV^ ■ Vpi = -VV*! 1 • Vp, 

A^i + e trace [D 2 i/j 1 D 2 (j)} + e 2 det D 2 ip! = pi - det D 2 <p. 

Differentiating the first equation with respect to space, we find the evolution 
equation for Vpi: 

^Vpi + ((V0 + eV^i) ± - V)Vpi 
= -(D 2 + eD 2 tp 1 )Vpi -D 2 i) l Vp 1 - -D 2 pV^i. (35) 
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We claim that in order to conclude the proof it is enough to have an estimate 
of the form 

11-01 OIIcmo?) < C(l + \\ Pl (t, -)llao.i (T2) ), (36) 

where C depends on 0. Let us admit this bound temporarily, and finish the 
proof of the theorem: using (}36|) and (J35|) . we obtain 

^||Vpi|| L ~ < C{t){\ + ||Vpi|| L - + e||Vpi|||oo), 

where the constant C(t) depends on the C 2 (T 2 ) norm of (p(t, •), 4>(t, •)). This 
quantity is bounded on every interval [0, T]. 

Thus we conclude using Gronwall's lemma that ||Vpi(i, OlU^Or 2 ) remains 
bounded on [0,T e ] with T € going to T as e goes to 0. We then choose T as 
large as we want, since when d = 2 the smooth solution to EI is global in 
time. From estimate (}36|) the W 1 ' 00 bound on pi implies a W 2 '°° bound on 
ip\. Then, we remember that 

P € -P yj, W ~ V0 

e e 
to conclude the proof of Theorem 16.31 

□ 



Proof of the estimate (I36j) 

We write the equation followed by ipi as follows: 

A^i = - trace [eD 2 ij^D 2 ^] - e 2 det D 2 ^ + pi - det £> 2 (/>. 
We recall that 

hence, using Schauder C 2 ' a estimates for solutions to Laplace equation (see 
[TB*]). we have 

Uxiy.a <d(l + e||^i||cv* +e 2 ||^i|lc^) ; (37) 

where Ci depends on ||0|| C 2,a, ||pi||c«- The inequality ([57)1 will be satisfied 
in two cases: either for ||^i|| C 2,a < C2 or for ||"0i I|o 2 =« > de -2 where 62,63 
are positive constants that depend on C\. 
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Now we show that ip e , solution of ()28|) . is bounded in C 2,a for p e bounded 
in C a norm. We consider for t G [0, 1] ip e t the unique up to a constant periodic 
solution of 

det(7 + e£> 2 ^) = 1 + tep e . 
Diiferentiating this equation with respect to t, we find 

MijDijdtipl = p e , 

where M is the co-matrix of I+eD 2 ^. From the regularity result of Theorem 
14. 3[ M is C a and striclty elliptic. From Schauder estimates, we have then 
\\dti>t \\c 2 ' a < C||p e ||c 2 . Q ; an d integrated over t G [0, 1], we get 

\m\c*-<c\\p e \\cp.«. 

Hence, since ip e = <f) + eipi, we have ipi bounded by C/e in C 2,a . Hence it can 
not be bigger than C 3 /e 2 , and to satsify (pT7J). we must have 

where C 2 as above depends on ||0|| C 2, a , ||pi||c a - This proves estimate (j3BJ). 

□ 
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